Electron-muon correlation as a new probe to strongly interacting quark-gluon plasma 
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As a new and clean probe to the strongly interacting quark-gluon plasma (sQGP), we propose 
an azimuthal correlation of an electron and a muon which originate from the semileptonic decay 
of charm and bottom quarks. By solving the Langevin equation for the heavy quarks under the 
hydrodynamic evolution of the hot plasma, we show that substantial quenching of the away-side 
peak in the electron-muon correlation can be seen if the sQGP drag force acting on heavy quarks 
is large enough as suggested from the gauge/gravity correspondence. The effect could be detected 
in high-energy heavy-ion collisions at the Relativistic Heavy Ion Collider and the Large Hadron 
Collider. 

PACS numbers: 25.75.Cj, 24.85.+p 



Quark-gluon plasma (QGP) is a novel state of matter 
expected to exist at extremely high temperature in early 
universe. Laboratory studies of the QGP by means of 
the high energy nuclear collisions are underway at the 
Relativistic Heavy Ion Collider (RHIC) and will be con- 
tinued at the Large Hadron Collider (LHC) 0]. The 
space-time evolution of the QGP created at RHIC is suc- 
cessfully described by ideal hydrodynamics where viscous 
effects arc assumed to be small [2J. This indicates that 
the quarks and gluons in the QGP are strongly interact- 
ing even above the critical temperature of the deconfinc- 
mcnt transition in accordance with previous theoretical 
expectations Q. The small value of a ratio of the shear 
viscosity rj to the entropy density s recently obtained 
from the gauge/gravity correspondence [i| and from the 
lattice quantum chromodynamics (QCD) simulations [H 
strengthens the idea of the strongly interacting QGP, or 
the sQGP in short. 

The hard probes such as high energy jets, heavy 
quarks, and quarkonia created in the initial stage of the 
heavy ion collisions would provide us with experimen- 
tal information on QGP. For example, suppression of 
the jet events and the disappearance of the jet corre- 
lations are intimately related to the parton energy loss 
inside QGP @, 0|- Suppression of the dileptons from 
the electromagnetic decay of the heavy quarkonia is re- 
lated to the screening of the heavy quark potentials 
in QGP Semileptonic decays of charm c and bot- 
tom b also carry information on the interaction of heavy 
quarks with QGP 0]: Recent data of the energetic sin- 
gle electrons from heavy quark decays 1^, 11 1 indicate a 
rather large en ergy loss of c and b inside the hot medium 
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The purpose of this Letter is to propose a new ob- 
servable which further signifies the dynamical property 
of heavy quarks in QGP: It is an azimuthal correlation 
of an electron (e) and a muon (fx) which originate from 
the semileptonic decays of back-to-back heavy quarks 
traversing the plasma. In particular, quenching of the 
away-side peak in the e-fx correlation turns out to be a 



clean observable to extract the in-medium energy-loss of 
heavy quarks [l^|. There are several advantages to look 
at the e-/x correlation over the other correlations: First 
of all, leptons from semileptonic decays keep the direc- 
tion of their parent heavy quarks as long as their trans- 
verse momenta are large enough (e.g. px > 3 GcV/c for 
electrons) [l9j]. Therefore, the dileptons reflect the corre- 
lation of a heavy quark-anti-quark pair just prior to its 
semileptonic decay. Furthermore, the e-fj, pair does not 
receive contamination from virtual photons and neutral 
vector mesons, so that the backgrounds are significantly 



reduced in comparison to the e + -e and 



correla- 



tions. The lepton-hadron azimuthal correlation [20), |21| 
could also provide us with useful information about the 
heavy quarks in the medium. However, the signal may be 
contaminated by other hadronic sources and interactions. 
Because of these reasons, the e-fi azimuthal correlation 
would be an ideal probe to study the heavy quark prop- 
agation in the sQGP. 

In the following, we describe the dynamics of heavy 
quarks by the relativistic Langevin equation with hy- 
drodynamic background, which is an approach recently 
developed by the present authors [l2j|. The hot matter 
created in the heavy ion collisions at RHIC is treated 
as a perfect fluid obeying the relativistic hydrodynamics 
equations 22]: 



0. 



T>™ = (e + P)u" 



Pg 
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where is the energy-momentum tensor, e is the local 
energy density, P is the local pressure, and w M is the 
local flow velocity. Initial conditions of the space-time 
evolution are determined to reproduce the experimental 
hadronic obscrvables in low px regions. 

The heavy quark dynamics in the hot matter is de- 
scribed by the relativistic Langevin equation formulated 
in the rest frame of a fluid element 13, 23|, 



At, Ap=-r(p)pAt + £(t), (3) 
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where x and p are the position and the momentum of a 
heavy quark with mass M, E(p) = ^Jp 2 + M 2 is the ki- 
netic energy of a heavy quark, T(p) is the drag coefficient, 
and is the Gaussian white noise with Dij(p) being 
the momentum dependent diffusion constant. Following 
our previous study, we adopt a particular parametriza- 
tion, T(p) = 7T 2 /M , where T is local temperature and 
7 is a drag parameter: This is motivated by the drag 
force obtained by using the gauge/gravity correspon- 
dence Then, the relativistic fluctuation-dissipation 
theorem leads to = 2~/T 3 (E + T)5 l3 /M 0. More 

sophisticated form with an anisotropic diffusion constant 
would be necessary in the future quantitative studies. 

Heavy quarks are produced through the initial hard 
processes simulated by the Monte Carlo event generator 
PYTHIA 6.4 [HI]. After the production, they are diffused 
inside the QGP fluid according to Eqs. (|3|) and (j4|) until 
the surrounding temperature decreases and reaches the 
critical temperature T c = 170 MeV. Below T c , heavy 
quarks hadronize to D or B mesons which subsequently 
undergo semileptonic decays simulated by PYTHIA. The 
color-singlet D and B mesons are assumed to propagate 
along the free-streaming path in the confined hadronic 
matter. 

Analysis of the nuclear modification factor Raa for 
energetic single electrons at RHIC on the basis of the 
Langevin + hydro approach leads to the drag parame- 
ter 7 = f-3 for heavy quarks 12j. This value is consis- 
tent with 7 ~ 2 estimated from the gauge/gravity cor- 
respondence applied to hot QCD matter [26J, while it 
is much larger than 7 ^ 0.2 given by the leading order 
perturbative QCD calculation [f3j|. If the heavy quarks 
were completely thermalized, the correlation of a heavy 
quark-anti-quark pair would be washed out. However, 
this is not the case as shown in our previous study [T3 | : 
Using the averaged temperature (T) ~ 0.21 GeV for the 
charm and bottom quarks in hot matter created by col- 
lisions at RHIC, one can estimate their relaxation time 
tq = Mg/-f(T) 2 . This together with the averaged stay- 
ing time (3-4 fm/c) of charm and bottom quarks inside 
the QGP fluid suggests that the charm quarks are par- 
tially thermalized while bottom quarks are not thermal- 
ized at RHIC for 7 = 1-3 12J . Therefore the correlation 



of a heavy quark-anti-quark pair is expected to survive 
partially during their propagation and serve as a good 
probe of the drag coefficient of heavy quarks. 

Let us now consider the e-fi azimuthal correlation 
obtained from our Langevin + hydro approach. The 
electrons and muons are from the semileptonic decays: 
D — > I, B — > I, and B — » D — > Z, where I stands for 
the lepton. Taking into account the current setups of the 
PHENIX detector, we take electrons with the pseudora- 
pidity \r)\ < 0.35 and the transverse momentum p^ > 3.0 
GeV/c as trigger particles, and muons with 1.4 < \rj\ < 
2.1 andp^ 1 ' > 3.0 GeV/c as associate particles. Collision 
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FIG. 1: (a) Azimuthal correlation of electrons and muons de- 
cayed from a heavy quark-anti-quark pair. Both the trigger 
electrons and the associate muons satisfy pr > 3.0 GeV/c. 
Trigger electrons are in mid-pseudorapidity < 0.35 and as- 
sociate muons are in forward-pseudorapidity 1.4 < \r)\ < 2.1. 
(b) Away-side quenching factor 7^ as a function of the drag 
parameter 7 with three cases of the transverse momentum cut 
{Pt — 0-3, 1-0, and 3.0 GeV/c) for associate muons. Frigger 
electrons satisfy ' > 3.0 GeV/c. Plots with the same drag 
parameter 7 are slightly shifted to avoid overlap. 



geometry is fixed to be nearly central with the impact pa- 
rameter b = 3.1 fm. 

Shown in Fig. [1] (a) is the e-fj, azimuthal correlation 
per trigger electron as a function of the relative open- 
ing azimuthal angle A0. If there is no energy loss of 
heavy quarks (7 — 0), high momentum muons appear 
in the away-side as a result of the back-to-back nature 
of the heavy quark production. It should be noted that 
the width of the peak in the 7 = case comes domi- 
nantly from finite primordial transverse momentum (~ 1 
GeV/c) of initial gluons inside a nucleon in the default 
parameter set of PYTHIA and that we do not consider 
nuclear broadening nor shadowing in the parton distribu- 
tion of a nucleus. There is only one peak in the away side, 
simply because a single semileptonic decay docs not pro- 
duce more than one lepton in the same direction. As the 
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energy loss increases, transverse momentum of the asso- 
ciate muons also decreases and the total yield in the away 
side with p^ > 3.0 GeV/c decreases. This is clearly seen 
for 7 = 3 and 10 in Fig. [I] (a). 

To quantify the effect of the suppression in the away- 
side muons and to compare the result with the future 
experiments, we introduce a quenching factor Iaa &s fol- 
lows: 



Iaa 



d(A(f>) 

^Aa/^PP: 



1 dN w 



iVtrig dA(j) 



ZYAM 



(5) 



where iVtrig and iV assoc are the number of trigger par- 
ticles (the electrons here) and the number of associate 
particles (the muons here). ZYAM implies zero-yield-at- 
minimum where muons other than the ones from heavy 
quark decays are subtracted as backgrounds. Theoreti- 
cally, Yiaa is a function of the drag parameter 7, while 
E pp (the case for pp collision) corresponds to 7 = 0. Tak- 
ing (0 mi n,0 max ) = (0,27r) in Fig. Q] (a), the away-side 
quenching factor I AA withp^ > 0.3, 1.0, and 3.0 GeV/c 
are evaluated as shown in Fig. [1] (b). We find that the 
quenching factor is sufficiently sensitive to 7: In the case 
of high Pt muons, I AA ^ s as small as 0.55 for relatively 
large drag parameter 7 = 3. We thus conclude that I AA 
is a useful observable which carries information of the 
heavy quarks in sQGP. 

Let us now consider the azimuthal correlation of an 
electron and a light charged hadron. The latter includes 
the charged pions, the charged kaons and the protons 
from the decay of the D and/or B-mesons. Although the 
detection of the electron-hadron (e-h) correlation itself is 
less demanding statistically, there are several disadvan- 
tages of the e-h correlation against the e-/i correlation: 
(i) light hadrons interact strongly with the hadronic mat- 
ter so that the azimuthal correlations are smeared by fi- 
nal state interactions, (ii) heavy quark propagation in 
the sQGP may affect the evolution of the flow of the hot 
matter and eventually introduces azimuthal angle depen- 
dence in the production of light hadrons [20j, |2l|, and 
(iii) for non-central collisions, anisotropic collective flow 
of the matter induces an e-h correlation not associate 
with the back-to-back correlation of heavy quarks. 

Assuming that the problems (i)-(iii) are not fatal, we 
make a theoretical estimate of the magnitude of the 
quenching in the e-h correlation. Here we adopt the 
same kinematic regions as those in preceding experiment 
[27l | in counting trigger and associate particles. We take 
electrons with transverse momentum 2.0 < p£ < 3.0 
GeV/c as trigger particles. Shown in Fig. [5] (a) is the 
e-h azimuthal correlation for the charged hadrons with 
0.5 < ' < 1.0 GeV/c. All the counted particles are in 
mid-pseudorapidity \r]\ < 0.35 and the impact parameter 
is taken to be b = 3.1 fm as before. In this case, there 
arise peaks in both near and away sides: The near-side 
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FIG. 2: (a) Azimuthal correlation of electrons and charged 
hadrons decayed from a heavy quark-anti-quark pair. The 
trigger electrons satisfy 2.0 < p^' < 3.0 GeV/c and the asso- 
ciate charged hadrons satisfy 0.5 < p^ 1 ' < 1.0 GeV/c. Both 
trigger and associate particles are in mid-pseudorapidity \r]\ < 
0.35. (b) Near-side and away-side quenching factors IXa as a 
function of the drag parameter 7 with two cases of the trans- 
verse momentum window (0.5 < < 1.0 GeV/c and 1.0 
f= Pt^ < 1-5 GeV/c) for associate particles. Trigger electrons 
satisfy 2.0 < < 3.0 GeV/c. Plots with the same drag 
parameter 7 are slightly shifted to avoid overlap. 



peak appears from the decay of a heavy meson into elec- 
tron and a charged hadron. As the heavy quark energy 
loss increases, both the near-side and away-side peaks are 
getting quenched. 

To quantify these quenching, we use the quenching 



factor Iaa as defined in Eq. 



with 



(0.57T, 1.57t) for the away side and (</> m i n , </>max) = 
(— 0.5-7T, 0.57r) for the near side. In Fig. [2] (b), near and 
away-side quenching factors I A ~ A with 0.5 < p^ < 1.0 
GeV/c and 1.0 < p^ < 1.5 GeV/c are shown. In the 
ideal situation where the problems (i)-(iii) can be ne- 
glected, the drag coefficient 7 = 3 leads to the near-side 
and away-side quenching factors as small as 0.35. It is 
thus not entirely hopeless to extract some information 
about the sQGP from I A ~ A although I AA allows us to 
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make much more transparent comparison between the 
theory and experiments. 

In summary, we have proposed an azimuthal electron- 
muon correlation as a new leptonic tool to study the 
energy loss of heavy quarks in the strongly interacting 
quark-gluon plasma. By using the Langevin equation 
of the heavy quarks under the background matter de- 
scribed by hydrodynamics, we have found that the away- 
side quenching factor I e AA is quite useful to quantify the 
energy loss of heavy quarks especially for leptons with 
high transverse momenta. We have also shown that the 
azimuthal clectron-hadron correlation could be valuable 
for studying the heavy quarks in hot matter, although 
there are various complications due to hadronic interac- 
tions. 

To correctly draw the information on the drag parame- 
ter inside the hot matter at RHIC, it is very important to 
determine the e-/j, azimuthal correlations experimentally 
not only for the nucleus- nucleus collisions but also for the 
proton-proton and proton-(deuteron-)nucleus collisions. 
Indeed, our quenching factor in Eq. ((5|) is defined as the 
ratio between the correlation in the nucleus-nucleus colli- 
sions and that in the proton-proton collisions. Moreover, 
studies on the width of the peak in the e-fj, azimuthal 
correlation in the proton-proton and proton-nucleus col- 
lisions may shed lights on the gluon distribution inside 
the energetic nuclei. The QGP at LHC will have higher 
temperature and longer life time than that at RHIC, so 
that the shape and magnitude of I AA would be modified 
even substantially. 

Finally, we mention possible theoretical improvements 
in the future: In our Langevin equation, we adopted the 
simplest form of the drag parameter and isotropic noise 
distribution. If the noise is not isotropic and depends on 
the direction of the heavy quark in the rest frame of the 
fluid as suggested by the gauge/gravity correspondence 
[24] ] , the shape and width of the azimuthal e-/i correlation 
could be affected. Possible interactions of the D and B 
mesons with the hadronic matter before the semileptonic 
decays, which is neglected in our calculation, may affect 
the magnitude of the away-side quenching. 
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